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Abstract 

We analyze the single-photon transport in a single-mode waveguide coupled to a whispering- 
gallery-type resonator interacting with a two-level atom. The single-photon transport properties 
such as the transmission and reflection amplitudes, as well as the resonator and the atom responses, 
are solved exactly via a real-space approach. The treatment includes the inter-mode backscattering 
between the two degenerate whispering gallery modes of the resonator, and the dissipations of the 
resonator and the atom. We also show that a generalized critical coupling condition, that the 
single-photon transmission at the output of the waveguide goes to zero on resonance for a matched 
system, holds for the full coupled waveguide-ring resonator-atom system. 

PACS numbers: 42.50.Ct, 42.79.Gn, 42.65.Wi 



* Electronic address: |jushen@stanford.edu| 
^Electronic address: lshanhui@stanford.edul 



1 



I. INTRODUCTION 



There has been a lot of recent interest in the systems of coupled whispering-gallery- 
type micro-resonators and atoms (a schematic configuration is 
shown in Fig. [T]). This configuration is particularly explored in a variety of fundamental 
and applied studies such as cavity quantum electrodynamics (cavity QED), single-atom 
manipulation and detection, and biochemical sensing. In the context of cavity QED, this 
kind of system has been treated semiclassically by assuming a weak classical input 10 1. 
Recently a comprehensive study of this system based on the semiclassical approach has 
been carried out by Srinivasan and Painter j^. Such semiclassical treatment is of direct 
current experimental interest, since the incident state in most cavity QED experiments are 
indeed weak classical beam. 

There are, however, several reasons for which the semiclassical treatment is not com- 
pletely satisfactory. First of all, given that these systems would ultimately be used to 
process quantum states of photons, it would seem valuable to treat these systems with fully 
quantized input states as well. Secondly, even in the presence of a classical input beam, 
the output of this system is not a coherent state. The semiclassical treatment typically 
provides various correlation functions of the output state. For a complete description of the 
output, however, one might instead prefer a direct description of the output quantum state 
itself. Thirdly, the semiclassical treatment is typically carried out numerically, by setting 
up a set of master equations for the evolution of the density matrix for the atom-cavity 
system, and by truncating the density matrix assuming a maximum number of photons in 
the cavity. For very weak classical input, for example, the maximum number of photon is 
often taken to be 1. For a deeper understanding of the underlying physics, however, exact 
analytic results are quite valuable. Finally, the physics intuition related to a semiclassical 
treatment has to be always treated with caution. For example, the truncation procedure, as 
discussed above, yields a set of coupled nonlinear ordinary differential equations, even when 
the maximum photon number in the cavity is fixed to be 1. Based upon this set of coupled 
nonlinear differential equations, one may be tempted to discuss the so-called "single-photon 
nonlinearity" . However, it is worthwhile to emphasize that nonlinearity has to arise from 
photon-photon interaction. 

Motivated by the discussions above, in this article, we consider the single-photon trans- 
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port for the system of coupled whispering-gallery-type micro-resonators and atoms, assuming 
a single-photon Fock state input. In contrast to the semiclassical treatments, where in most 
cases only numerical results are obtained so far, here the use of a fully quantized formalism 
allows one to straightforwardly obtain analytic and exact results of the transport properties 
and the system responses, which has not been done before for this system. Our analytic 
results are in well agreement with the experimental data js, [ll| and with the numerical 
results in Srinivasan and Painter's for the parameters considered in their studies j^. Us- 
ing these analytic results, we systematically explore the parameter regime that is of direct 
experimental interests, to understand the delicate interplay between the various types of 
coupling and the intrinsic atom and resonator dissipations. 

This article is organized as follows. In Sec. [Illwe first introduce the system configuration 
and the Hamiltonian. In Sec. IIIIl the single-photon transport properties and the system 
responses are then solved analytically. Sec. [IV] describes the general properties of the exact 
solutions, such as spectral symmetry of the transmission spectrum, and generalized critical 
coupling condition of the full coupled system. In Sec. |Vl we compare our results with some 
experimental results of waveguide-resonator systems, and with published numerical results 
of coupled waveguide-resonator-quantum dot systems. Sec. |Vl] presents a systematic study 
of the single-photon transport of the full coupled system of waveguide-ring resonator-atom. 
Sec. IVIII discusses the effects of the intrinsic dissipations of the resonator and the atom. In 
Sec. IVIIII we investigate the effects of complex inter-mode backscattering. Sec. HXl discusses 
the effects of detuning between the resonator and the atom. Finally, Sec. |X] sums up the 
article, pointing out some applications and generalizations of the formalism. 



II. THE SYSTEM AND HAMILTONIAN 



The system of interest in this article is schematically shown in Fig. ([T]): a whispering- 
gallery type resonator interacting with a two-level atom is side-coupled to a single-mode 
waveguide. A whispering-gallery type microresonator, such as a ring resonator, a micro- 
sphere, a microtoroid or microdisk, supports two degenerate whispering gallery modes 
(WGMs) that propagate around the resonator in opposite directions. (Throughout this 
paper, we will often use the term "ring resonator" to designate a whispering-gallery type 
microresonator). We also include the interactions of the atom and the resonator with the 
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reservoirs. Such interactions with reservoirs give rise to intrinsic dissipations [id. Il2[|. 
The Hamiltonian of the composite system S @ R is H = Hs + Hr + Hsr- 

+ j dx5{x)[VaC^R{x)a + V*a^CR{x)] + j dx 5{x)\Vbc\{x)h + V^l)^ cl{x)] 

+ [hb^^a + /i*a"f6) , (la) 

HR/h=J2 ^lAfij + ^^1^1 f 2^ + ^3ir5/3j , (lb) 
j j j 

HsR/h = J2 {^Hj^ + ^j^^^h) + {^Hj^ + s-^^^2j) + J2 {^Hj"^- + ^j^+^3j) • 
j j j 

(Ic) 

if5 is the Hamiltonian of the system S of coupled waveguide-resonator-atom. This Hamil- 
tonian includes the waveguide, the resonator, and the atomic part, as well as the interaction 
between the waveguide and the resonator, and the atom and the resonator. 

The first line in Hs of Eq. ([Taj) describes the propagating photon modes in the waveguide. 
c^y^(x) is a bosonic operator creating a right /left-moving photon at x. is a reference 
frequency, around which the waveguide dispersion relation is linearized [3]. Its value does 
not affect the photon transport properties. 

The second line in Hs describes the modes in the resonators and the atomic states. For 
the resonator, is the creation operator for the counter-clockwise WGM mode and is 
the creation operator for the clockwise WGM mode, both of frequency Uc- For the atom, 
a^g{al) is the creation operator of the ground (excited) state, cr+ = aja(,(cr_ = ajoe) is the 
atomic raising (lowering) ladder operator satisfying a+\n,nc = 0,—) = \n,nc = 0,+) and 
cr+|n. Tie, +) = 0, where |n, ric, ±) = \n) (g) \nc) ® |±) describes the state of the system with n 
propagating photons, photons in cavity mode, and the atom in the excited (+) or ground 
(— ) state. Qe — ^g{= ^) is the atomic transition frequency. 

The third line in Hs describes the interactions between the waveguiding modes and the 
WGMs. Va/b is the waveguide-resonance coupling strength of each WGM. Here the right- 
moving (left-moving) mode only couples to the phase-matched counter-clockwise(clockwise) 
mode. 
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The fourth hne in Hs describes the interactions between the atom and the whispering 
gallery modes, ga and gb are the resonator-atom coupling strength for each respective WGM. 

The last line in Hs describes the inter-mode backscattering between the two degenerate 
WGMs, induced through imperfection of the resonator, h is the inter-mode backscattering 
strength. 

Hr of Eq. (jlbp describes the reservoir, which is composed of three subsystems: R = 
-Ri©-R2©-R3- -Ri, R2 and -R3 are assumed to be independent. Each _Ri, R2 and -R3 
is modeled as a collection of harmonic oscillators with frequencies uJij, u!2j, and u^j, and 
with the corresponding creation (annihilation) operators r\j (rij), rlj (r2j), and rlj {r^j), 
respectively. 

HsR of Eq. fITcl) describes the interactions between the resonator and the atom with the 
reservoirs, respectively. The WGM couples to the jth reservoir oscillator rij in Ri, while 
the WGM couples to the jth reservoir oscillator in i?2. Both WGMs are assumed to 
couple with the reservoir with the same coupling constant kj. The atom cr+ couples to the 
jth reservoir oscillator r^j in R3 with a coupling constant rjj. 

By incorporating the excitation amplitudes of the reservoir R, it can be shown that the 



effective Hamiltonian H^s of S can be obtained and is given by 



13|: 



+ J dx 5{x)[VaC^R{x)a + V:a^CRix)]+ j dxS{x)[Vbc[ix)b + V,*b^CL{x)] 
+ {9aaa+ + fi'>V_) + {gbba+ + glh^crJ) 

+ {hh^a + h*a%) , (2) 

where 1/tc = 7c and 1 /r^ = 7^ are the intrinsic dissipation rates of the resonator WGM 
and the atom, respectively, due to coupling to the reservoir. We will call iiTefr as H in the 
following for brevity. Note that \Va/h\'^ /vg, 1/tc, l/r^, ga, gb, and h all have the same unit 
as frequency. 

Here we make some remarks on the Hamiltonian of Eq. ([2]) : 

1. The coupling constants are determined by the underlying photonic and electronic 
states of each constituent of the system. The relations between these coupling con- 



stants are not arbitrary but rather are constrained by symmetry of the system. These 
symmetry properties in turn are reflected in the response spectra of the system. The 
relevant symmetry operations here are mirror and time-reversal symmetries. A dis- 
cussion of the symmetry transformations is given in Appendix [XI 

2. Except obeying the symmetry transformations, the specific numerical values of the 
coupling constants depend upon each specific experimental realization. A detailed 
discussion is provided in Appendix [Bl Here, we just note that that in general \ga\ 7^ \gb\, 
and both Qa and gi, can be complex numbers. 

In the following, we first derive the exact solutions to the single-photon transport and 
the system responses {i.e., WGMs and atom excitations) described by the Hamiltonian of 
Eq. ([2]), without placing any constraints on the coupling constants. Later, we will specialize 
to specific choices of the coupling constants to compare with experimental data, and to some 
published numerical results. 

III. EXACT SOLUTIONS FOR SINGLE-PHOTON TRANSPORT 

We consider the temporal evolution of an arbitrary single- photon state as described 

by the Schrodinger equation 



where |0) is the vacuum, which has zero photon and has the atom in the ground state. 
(f)R/L{x,t) is the single-photon wave function in the R/L mode. ea/b{t) is the excitation 
amplitude of the whispering gallery mode, and eq{t) is the excitation amplitude of the atom. 
For this state, the Schrodinger equation (Eq. ([3])) thus gives the following set of equations 




(3) 




+ e,(t)at|0) + e,(t)6t|0) + g^(t)at|0) 



(4) 
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of motion: 



d ~ „ ~ d ~ 

-iVg—(t)R{x, t) + 5{x)Vaea{t) + (u;o + fig) <Pr{x, t) = l—(j)n{x, t), 

d ~ ^ ~ d ~ 

+iVg—(f)L{x, t) + 6{x)Vbeb(t) + {uo + Qg) (I)l{x, t) = t), 

1 ~ d 

{u, + ng- i-)ea{t) + KVr(0, t) + glB.it) + h*eh{t) = «^e,(t), 

1 „ ~ „ . (9 „ 

{uc + ng- i-)eb{t) + V*(f)L{0, t) + gleg{t) + = ^^eb(t), 



1 . . _ d ^ 

{Qe - ^ — )(iq(t) + gaCait) + 5'fee6(t) = 2 — eq(t). 
To CI 



(5a) 
(5b) 
(5c) 
(5d) 
(5e) 



For any given initial state |$(t = 0)), the dynamics of the system can be obtained directly by 
integrating this set of equations (Eqs. (JSj)). In this way, one could study the time-dependent 
transport of an arbitrary single-photon wave packet. 

In the following, we concentrate on the steady state properties. When |$(t)) is an eigen- 
state of frequency e, i.e., |$(t)) = e~*''*|e"'"), Eq. ([3]) yields the time-independent eigen 
equation 

H\e+) = he\e+). (6) 

and the interacting steady-state solution \e~^) can be solved for. Here he is the total energy 
of the system. 

For an input state of one-photon Fock state, the most general time-independent interact- 
ing eigenstate for the Hamiltonian of Eq. ([2]) is: 



dx 



</.r(x)4(x) + (f)L{^)c[{x)] 10) + e,at|0) + e,6t|0) + e,a+|0), (7) 



where we denote the time-independent amplitudes by the corresponding untilded symbols, 
e.g. ea{t) = Ca e"*"^*, etc. The connection between the interacting eigenstate and a scattering 
experiment is described by the Lippmann-Schwinger formalism jl4j. Il5l. Il6[|. 

The time- independent Schrodinger equation of Eq. (P) for the state |e+) of Eq. yields 
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the following equations of motion: 

d 

-Wg—(f)R{x) + 5{x)Vaea = {e-uo- fig) (j)Rix), (8a) 
d 

+Wg—(t)L{x) + 5{x)Vheh = {e-ujQ- Vtg) (piix), (8b) 
(^c - «-)ea + V*4>r{0) + gleg + h*eb = (e - Qg) Ca, (8c) 
{uc - i—)eb + Vj*0l(O) + gleg + hea = (e - fig) e^, (8d) 
(fi - i—)eq + ^faCa + ^'foCfe = (c - fig) 6^, (8e) 

with e = + fig, and u = ujo + VgkR. Our aim is to solve for the transmission and reflection 
amplitudes for an incident photon. For this purpose, we take 4>r{x) = e**^^' {6{—x) + t6{x)), 
and 0l(x) = re~'^^^6{—x), where t is the transmission amplitude, and r is the reflection 
amplitude Q- The set of equations of motion, Eqs. fl8a|) -f l8e|) now read: 

-iVgit-l) + Vaea = 0, (9a) 
+Wg(-r) + Hefe = 0, (9b) 

{uc - i—)ea + K*-^ + 9*a(ig + h*eb = uiea, (9c) 

1 T 

{uc - )efe + V;*- + gleg + hea = ^e?,, (9d) 

Tc 2 

{Q-i—)eg + gaCa + gbGb = i^Cg, (9e) 
which can be solved straightforwardly for Q, t, r, e^ and e^: 



U — Uq 



(10a) 



[u~u)c+ ^(uj -n + i:^'j (^LLi-LLic+iy-'j - G^] + (tJ - H + i:;^) - g*gi,/i - g^g*/!* - [/ip (a; - Q + i^^) + iG^ T 

(uj-ujc + i^+ ir) [(w - n + i-i) (a; - a;c + i-;^ + - G2_] ~ gjgj/i - gag^h* ~ \h\^ (uj - Q + i-i) 

(10b) 

-i^^ [gag;+h(u-n + ^^)] 





-n + 




- cJc + — 


+ ir 




- 9l9bh - 


-9ag^h* - 






^9 J 










— Ll>c - 


\-i^+iV^ ~ \gb\- 


'] 








(uj - ujc + i:^ + iT^ [(<^ - 


- n + 






+ iT 




- 9a9bh - 


'9a9th* ' 




(a; - n H 








a 


[dagt + h 






)] 












- n + 




- 0Jc + i-f- 


+ ir 




-9*a9bh~ 


-9a9^h* ' 




(w -Q-\ 


^qJ 








(o; — uic - 




+ ir) + 













(lOc) 



(lOd) 



^6 = 7 ^ TVTT : T^-N -TT^ ^ —77 1 TX' (^Oe) 



(uj-0Jc+i:^+ ir) [(oj - n + j;^) (w - + 4:;]^ + - G^ j - gjg^/i - gag^h* - \h\^ (u; - Q + i^) ' ^ ^ 
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where we have used IVJ = \Vb\ =V,T = G\ = \ga\ +\gb\ , and = lotl — Iflal , all are 
real numbers. Notice that F is the external linewidth of the WGM's due to waveguide-cavity 
coupling. These analytic expressions of the amplitudes, Eq. (]10bp - (]10f|) . provide a complete 
description on the single-photon transport properties. These equations are applicable to 
arbitrary two- level systems having any orientation of the electric dipole moment, specified 
through Qa and g^. 

IV. GENERAL PROPERTIES OF THE TRANSMISSION SPECTRUM 

The expressions of the amplitudes (Eqs. ffTOj) ) allow us to make exact statements on 
the general properties of the spectrum of the full coupled system, such as the spectral 
symmetry properties and the generalized critical coupling conditions. These statements are 
valid even in the presence of resonator and atom dissipations. We will exemplify these 
general properties with concrete examples in the following sections. 

A. Spectral Symmetry Properties of The Transmission Spectrum 

It is straightforward to show that the transmission amplitude of Eq. (llObI) satisfies the 
following spectral symmetry condition, regardless of the values of the cavity and atom dis- 
sipations: 

t* {-6uj,-A,-6e) =t{6uj,A,6e), (11) 

where 6uj = u ~ Uc, A = Q — Uc is the frequency detuning between the resonator and the 
atom, and 66 = 6h + 0^ — 6a — is the phase mismatch. This condition states that the 
transmission spectrum T(u;) = of one system with parameters (—A, —66) is mirror- 

imaged with respect to Uc to that of another system with parameters (A, 66). 

In particular, for a system with A = and 66 = (mod 2tt), the transmission spectrum 
is symmetric with respect to u = Uc'- T{uj = ujc + 6uj) = T{uj = ujc — 6uj). The former 
condition A = says the resonator and the atom are in-tuned {Vt = ujc); while the latter 
condition 66 = (mod 27r) requires g^gth to be purely imaginary so that gagbh + gaglh* = 0. 
The same spectral symmetry property holds for other amplitudes of Eq. f llOcp -( fTOfj) . 
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B. Generalized Critical Coupling Condition 



Critical coupling is defined for a system of coupled waveguides and resonators when 
transmission of the input signal goes to zero at the output port at resonance. The existence 
of critical coupling is of direct interest to quantum optics experiments where one often desires 
to eliminate single-photon transmission so that the signatures of two-photon transmission 
are distinct from single-photon effects. Here we present the condition under which critical 
coupling can be achieved for single-photon transport in the full coupled waveguide-ring 
resonator-atom system. 

By requiring the transmission amplitude of Eq. (llObp to be zero, one can show that the 



critical coupling is reached aX u = Uc if the following criterions are satisfied: 
1. The resonator and the atom are in-tuned, i.e., Q — Uc = A = 0; and 
2. 

r' + r,G'r=( l-Y + TlGl+\h\', (12a) 
glQhh + gaglh* = Q. (12b) 



Eq. fll2ap is a magnitude-matching condition. Eq. (112bp is a phase-matching condition, 



requiring g^Qhh to be purely imaginary. Since gl^Qhh = i\ga\\gb\\h\e^^^ , with 69 = 9^ + 
9b — 9a — 71/2, Eq. (112b[) is equivalent to the phase-matching condition 69 = (mod 27r). 

V. NUMERICAL VALIDATION 

Before we proceed to understand the detailed predictions of Eqs. ( fTOl) . here we first 
compare Eq. (110bl) -( fT0fj) to experimental data of waveguide-resonator systems j^, 11], and 



to numerical results of Srinivasan and Painter's [6|] on the waveguide-microdisk-quantum dot 
system. In all cases the agreements are excellent. 

A. Fitting to Experimental Data of Waveguide-Resonator Systems 

To demonstrate the validity of the exact solutions, we apply these expressions to fit two 
of the recent experimental data of waveguide-microtoroidal resonator systems 

The 
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results are plotted in Fig. ([2]). The amplitudes for the waveguide- resonator system can be 
deduced from those of the full coupled waveguide case by decoupling the atom. In particular, 
the transmission amplitude is given by 

tM = ) , . (13) 

The details are given in Appendix O Note only the magnitude of h appears in the amplitude, 
its phase does not. 

F,g. 0(a) show. ..a,.,„.s,o,> spect.a of a wavegu,de-.ic.oto..oida. .esoua^o. .yste,„ g 

In the experiments, the waveguide-resonator coupling strength |^| is varied by adjusting the 
distance between the waveguide and the microtoroid. The lower trace is taken such that the 
critical coupling condition is approximately satisfied: ^ \hf + (I/tc)^; while the upper 
trace is for conditions of under-coupling: F^ ^ 1^1^ + (V'^'c)^- By fitting to the spectra using 
Eq. (IT^ . one can extract the numerical values of \h\, 1/tc, and F. Our numerical results 
show that \h\ ~ 1/tc. 

When the same values of \h\ and I/tc determined from fitting to the lower trace are used 
to fit the upper curve, the width of the spectrum is narrower than that of the experimental 
data, and the minimum region has significant deviation. If, however, somewhat larger 
values of \h\ and 1/tc are used instead to fit the data, the spectrum can be fitted reasonably 
well, as shown by the upper blue curve in Fig. ([2]) (a). This indicates that both the inter- 
mode backscattering and the intrinsic loss of the resonator are slightly suppressed when the 
coupling between the waveguide and the resonator is strong. 

As another validation, in Fig. ©(b), we use Eg. ffT3l) to fit the transmission spectrum of 
another waveguide-microtoroidal resonator system 11[, which is in the deep under-coupling 
regime such that F^ <^ \hf + (I/tc)^. The numerical fitting indicates that \h\ ^ 1/tc ^ F 
so the inter-mode backscattering dominates. In this regime, the transmission spectrum 
shows a doublet structure. The doublet structure in the transmission spectrum induced 
by the inter-mode backscattering for ultra-high-Q resonator in the > F regime is well- 
known jll . 



19 



20|. 



In both cases above, the extracted numerical values agree reasonably well with those 
reported in the experiments. 
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B. Comparing to Numerical Results of Waveguide-Resonator- Atom Systems 



We next compare our analytic results of the amplitudes of Eqs. (110b|) - (]10f[) to the 



numerical results on coupled waveguide-microdisk-quantum dot system in Srinivasan and 
Painter [6|, where gi, = Qa by properly choosing the orientation of the quantum dot dipole 
polarization, and the azimuthal origin. From Eqs. (fT0bl)-(fT0l. one immediately sees that 
all amplitudes, except e^, depend on the magnitude of Qa only but not its phase. The atom 
excitation amplitude Cq is proportional to the spectrum \eq\ however also depends on 
\ga\ only. 

Fig. [3] plots the transmission and reflection spectra using the parameters corresponding 



to those in Fig. 5 of Ref. 



6|. (The relation between the parameters used in this article and 



Ref. is: 5f = l^fal = | -igo\, T = Ke, l/r^ = Ki = kt - ^e, l/r^ = 7|| +7p, and h = -13). In 
Fig. [3l the upper panel plots the spectra when the quantum-dot dephasing 7^ is zero; while 
the lower panel plots the spectra when 7^ ^ 0. Each panel plots the in-tuned [Vt = Uc) and 
de-tuned {Q 7^ ujc) cases. The results of both approaches are in excellent agreement. We 
were able to obtain excellent agreements with all results in Ref. 6| (not shown here). 

In Fig. [3] we also plot the atom excitation. In general, there are only two resonances 
in the atom excitation, in contrast to the three resonances in the transmission spectrum. 
Since the atom is an intermediary for correlating photons, one expects the photon-photon 
correlations are qualitatively different at the three resonances of the transmission spectrum. 
This is supported by comparing the atomic excitation spectra with the photon correlation 
function g^'^^r) in Fig. 12 of Ref. ^. 

VI. SINGLE-PHOTON TRANSPORT OF THE COUPLED WAVEGUIDE- 
RESONATOR- ATOM SYSTEMS 

We now turn the attentions to the single-photon transport of the full coupled waveguide- 
resonator-atom system and the system responses. The analytic results enable us to present 
a systematic parametric study of the transport properties by continuously varying the pa- 
rameters. Such a systematic study facilitates the understanding of the underlying physics 
in this complicated coupled systems. As a concrete example, we will take Qa = gb, and let 
g = \ga\ throughout the discussions hereafter. This choice of coupling constants correspond 
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to that considered in Ref . (6|] . Other choices of the couphng constants could be investigated 
in the same manner. 

In this section the inter-mode backscattering strength h is taken as real and the system 
is assumed lossless (1/tc = l/r^ = 0). The atom and the resonator are also assumed to be 
in-tune (f2 = u^. The effect of dissipation, complex /i, and atom-cavity detuning will be 
considered in later sections. 

We consider the spectra of transmission T = (Fig. H]), the group delay d(j)/du} (Fig. [5]), 
the atom excitation |eg|^ (Fig- E]), and the phase matched WGM excitation |ea|^ (Fig. [7]) 
for different value of g and h. (For completeness, we also plot the spectra of reflection 
R = |r(ci;)p and the counter-propagating WGM excitation le^p in Appendix [D|) . To facilitate 
visualization, these spectra are presented in a matrix form with different rows or columns 
corresponding to different g and h, respectively. Comparing this set of spectra provides 
very useful information to understand the full coupled system, since the atom and the 
WGM excitation determine the nature of the transmission resonances. Below, we refer to a 
resonance with large (small) atom excitation as a resonance of atom (cavity) nature. 

A. First column: the atom decoupled 

The leftmost column describes the waveguide-resonator case with the atom decoupled: 

1. When h = 0, the resonator acts as an all-pass filter with transmission T = 1 for all 
frequency. The group delay shows a Lorentzian peak. 

2. When h is slightly increased but small than F, a transmission dip develops at the 
resonance frequency uj = Uc, with 



Also, the transmission has a non-Lorentzian lineshape. 

3. When h = T, the transmission spectrum has a fiat bottom centered at uj = ujc where 
the transmission becomes zero. The transmission has a maximally-flat 2nd-order But- 
terworth filter lineshape given by 




(14) 




(u; - Uc) 



(15) 



{u-Uc)^ + 4\h\ 



4' 
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with a full width at half maximum (FWHM) equal to h \ . The group delay however 
shows two sphtted peaks. 

4. When h >T, the transmission shows two resonance dips at which the transmission is 
zero. The spectral separation between the two dips is 2-^|/ip — F^, which approaches 
2\h\ when \h\ > T. 

Thus the qualitative behavior of the transmission spectrum is determined by the ratio of 
|/i|/r. For an ultra-high Q resonator with very weak external coupling to the waveguide, 
i.e., very small value of F, even a weak mode cross-talking would qualitatively change the 
transmission spectrum. 

Since the atom is decoupled, all the transmission dips (resonances) have a cavity nature, 
i.e., the atom excitation is zero. 

B. First row: no inter-mode backscattering in the resonator 

The uppermost row describes an ideal resonator without mode cross-talking, i.e., h = 0: 

1. When (7 = 0, the atom is decoupled, and the resonator is an all-pass filter as discussed 
above. 

2. When g is slightly increased from zero such that 2g'^ <^ F^ (shown in figure: g = 0.3F), 
a transmission dip down to zero occurs at frequency u = uJc. This transmission dip is 
induced by the atomic resonance and is in contrast to that induced by small h above, 
which does not go to zero. In this regime, the poles of the transmission amplitude 
are u — ujc = —i^r, — ^), and —iT. The transmission spectrum can be well 
approximated by a Lorentzian with FWHM equal to Ag'^/T. Small deviation of the 
spectrum from a Lorentzian occurs only when \u! — u^d — F, where the poles with 
larger imaginary parts become important in the spectral response. Moreover, the 
atom excitation has a strong peak at u; = oOc, thus the transmission dip has an atomic 
nature. The group delay shows a single peak. Notice that the group delay in this case 
is much larger than case Bl above. 

3. When g is further increased such that 2g^ = F^, the transmission spectrum has a flat 
bottom centered at u = Uc where the transmission is zero. The transmission has a 
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maximally-flat 3rd-order Butterworth filter lineshape given by 

T{uj) = , (16) 

(u; - ujcr + Sis'! 

with a full width at half maximum (FWHM) equal to 2V2g. Both the group delay 
and the atom excitation however shows two splitted peaks. 

4. When g is further increased such that 2g'^ > (shown in figure: g = 2.5r), three 
transmission dips develop. This is because the two counter-propagating WGMs of the 
resonator, through linear superposition, can alternatively be described as two standing- 
wave modes. One of the standing-wave modes has zero amplitude at the atom location, 
and corresponds to the middle transmission dip at = cj^. This resonance is of cavity 
nature. The other mode has non-zero amplitude at the atom location, and experiences 
Rabi-splitting through interacting with the atom. The transmission spectra exhibits 
three dips when such Rabi splitting is large enough compared with F. The two side 
dips, which result from the Rabi splitting, have a mixture of cavity and atom nature, 
as confirmed by examining the atom excitation plot. The spectral separation between 
the two side dips is 2^2\g\'^ — F^. In addition, the two side-dips have a larger group 
delay than that of the middle peak, since the photon experiences delay from both the 
cavity and the atom. 

5. When 2g^ ^ F^ (shown in figure: g = 8F), the spectral separation between the two 
dips approaches 2\^\g\. The maximum transmission between each side-dip and the 
middle dip approaches 1 as 1 — XtTI^- Moreover, the width of the middle dip is the 
sum of the widths of the two side-dips. 

The case of h = represents the ideal case for WGM resonators and therefore is of 
fundamental importance. In Appendix [E], we summarize all analytic results related to this 
important case. 



C. General case: non-zero h and g 

We now discuss the general features of the cases when both h and g are non-zero. 

1. Comparing the leftmost two columns in Fig. HI where we increase g slightly from 
zero to a small non-zero value (0.3F in the figures). For all values of h, the effect of 
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introducing an atomic resonance by a small g is to create a resonance of atomic nature 



at = uJc, resulting in an asymmetrical Fano-type lineshape in transmission [2l|, 
and a large group delay at resonance. The rest of the spectrum is not significantly 
perturbed. 

2. Comparing the uppermost two rows in Fig. HJ where we increase h slightly from zero 
to a small non-zero value (O.SF in the figures). For all values of (? 7^ 0, the effect of 
introducing an inter-mode backscattering by h is to slightly distort the transmission 
spectrum to be asymmetrical. 

3. The cases with intermediate values of g and h can be obtained from perturbing from 
its counterpart in each direction. The exact lineshapes have to be computed using the 
analytic expressions. 

4. One important feature is the anti-crossing between the atomic and cavity resonances 
as one varies the values oi h 01 g. We will demonstrate this point using the column of 
g = 2.5F. Fig. [8] plots the real part of the poles of the transmission amplitude 
which indicate the spectral locations of the resonances. At /i = 0, the transmission 
spectrum starts with three resonances with the middle dip of cavity nature, and the 
two side-dips mixture of both cavity and atom nature, as previously explained. When 
h is increased (from h = O.SF to 10 F), the left two dips anticross. During the anti- 
crossing process, the nature of the resonance of the two dips is exchanged, as can be 
seen by examining the weight of the atomic excitation along the evolution process. In 
the end, when h becomes lOF, the left resonance is of purely cavity nature, while the 
middle resonance is of purely atomic nature, and the right resonance is of largely cavity 
nature with a little atomic mixture. Similar anti-crossing behavior can be observed 
when g is continuously varied instead. Thus, it is difficult to ascertain the nature 
of a resonance judging from the position of the resonance alone, in the case where 
significant back scattering occurs. 



VII. EFFECTS OF DISSIPATIONS 



The effects of intrinsic dissipations on the transmission resonances strongly depend upon 
the nature of the resonances. When a resonance is of purely cavity nature, the transmission 
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at the resonance frequency is insensitive to the atom dissipation, but only to the cavity 
dissipation. Similarly, when a resonance is of purely atomic nature, the transmission at the 
resonance frequency is insensitive to the cavity dissipation, but only to the atom dissipation. 
For a resonance of mixed nature, the transmission is affected by either type of dissipation. 

This is clearly seen in all cases in Fig. [HI where we introduce dissipations to all cases 
considered in Fig. HJ which are lossless. We choose dissipations that are either purely from 
the atoms (l/rg = F, I/tc = 0), or from the cavities (I/tc = F, l/r^ = 0). 

As a more detailed example. Fig. [10] plots the effects of increasing dissipations for the 
same case oi h = lOF and g = 2.5F. The left resonance, which is of purely cavity nature, is 
essentially unaffected by even very large atom dissipation. The right resonance, which has 
a little atomic nature, has its transmission minimum gradually lifted and its width slightly 
broadened at large atom dissipation. The cavity dissipation, on the other hand, strongly 
affects all three resonances, since the phase matched WGM has weights at each resonance. 
Among the three resonances, the middle one is least affected, since it is primarily atomic. 
At large cavity dissipations, it still leaves a small signature in the transmission spectra while 
the contributions from the other two resonances are no longer visible. Thus, the nature of 
the resonance influences the properties of the transmission spectrum even in the regime of 
large dissipation. 

VIII. THE EFFECTS OF COMPLEX h 

When Qa = Qb, all the amplitudes in Eq. (110bl) - (110fj) depend upon the phase of h only but 
not of Qa- In the above discussion, for concreteness, we have assumed h to be real. Here we 
consider the effects of a complex h on the transport properties. 

With Qa and Qb real, the spectral symmetry properties mentioned in Eq. ffTTj) in Sec. IIVI 
now reads 



Fig. [TT] shows a series of transmission spectrum with g = 2.5F, and h = 2.5e*^''F and 
A = 0, i.e., the cavity and the atom are in-tuned: 

1. When 6h = 0, two resonances of very different quality factor Q are close to each other 




(17) 



where 66 = 6, 



^. The spectral symmetry is determined by the phase of h only. 



at u; < ujc- 
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2. As 9h = 7r/4, the spectrum has three dips and is asymmetric with respect io uo = uj^. 

3. When = 7r/2, satisfying the condition that 6h — '26 g = f + nir with n an integer, 
the spectrum is symmetric with respect to u = Uc- 

4. The spectrum of 6^ = Sir /A is the mirror image of the spectrum of 6^ = tc /A with 
respect to u = Uc- 

5. Finally, the spectrum of 6h = tt is the mirror image of that of 6h = with respect to 

UJ = LUc- 

A general case is plotted in Fig. [121 with atom-cavity detuning A 7^ 0, and with finite 
atom and resonator dissipations. The transmission spectra with (—A, —6) and with (A, 6) 
are mirror-imaged with respect to u = Uc- 

IX. EFFECTS OF DETUNING {n / tOc) 

In this section, we discuss the effects of detuning between the resonator and the atom, 
with h = 0. Fig. [13] plots the transmission spectrum for the lossless case as (7 or A is varied. 
We note that the spectrum is always asymmetrical with respect to u = Uc when both A and 
g are non-zero. We now discuss these spectra plots, as organized in a matrix form. 

A. The first row: large detuning with |A| » F 

1. When (yf ^ A ((yf = 0.3F as plotted), the atom is decoupled from the system and 
creates a narrow atomic resonance aX u = Q. The background is the transmission of 
waveguide-ring resonator subsystem that has unity transmission. 

2. When g is increased {g = \j\f2Y as plotted), the atom becomes weakly coupled to 
the waveguide-ring resonator subsystem. Analytically, such a weak coupling regime 
occurs in the parameter range 2\/2g < A, where 2\/2g is the Rabi frequency, and 
\/2g < F. In this weak coupling regime, the transmission spectrum exhibits an atomic 
resonance dip aX uj = Vt. Also, the weak coupling of the atom and the cavity results in 
a small scattering between the two WGM, and consequently a small dip is present at 
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the resonant frequency uj = ujc- Notice that the transmission dip at the cavity resonant 
frequency does not reaches zero at its minimum. 

3. For large g such that 2^/2g > A and ^/2g > T, the atom is strongly coupled to the 
waveguide-resonator subsystem, and the transmission spectrum develops into three 
resonance dips, each reaching zero at its minimum. . 

B. Second column: decoupled and weakly coupled atom 

For decoupled {g <^ A) and weakly coupled {2\^g < A, ^/2g < F) atom, the atom 
creates a narrow atomic resonance at u = Q, and the location of the resonance moves with 
A. At A = 0, the weakly coupled atom mixes with the cavity resonance so the line width is 
slightly increased, compared with that of the A 7^ cases. 

For intermediate values of g and A, the exact lineshape and the weight of the nature 
of a resonance has to be computed using the exact expressions. Fig. [H] plots the effects 
of intrinsic losses of the resonator and the atom to the transmission spectrum. Again, the 
resonator loss strongly suppresses a resonance of cavity nature, and the atomic loss strongly 
suppresses a resonance of atom nature. 

X. SUMMARY 

We have provided a full quantum mechanical approach to treat the coupled waveguide- 
ring resonator-atom system, and derived the analytic solutions for the single-photon trans- 
port. The real-space approach outlined in this article can be generalized straightforwardly 
to treat cases such as cascaded multi-ring resonator, multi-atom, or multi-port configuration 
that is relevant to applications of add-drop filter, single-photon switching and delay lines. 
Our formalism can also provide a starting point for treatment of pulse propagation in these 
more complicated systems. 
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APPENDIX A: SYMMETRY TRANSFORMATIONS OF THE HAMILTONIAN 

The relations between the couphng constants are not arbitrary but rather are constrained 
by the symmetry of the Hamiltonian. The symmetries of interest for our configurations 
are mirror symmetry and the time-reversal symmetry. In this section, we consider these 
symmetries that constrain the form of the Hamiltonian, and establish the relations between 
the coupling constants. Note that the specific form of the relations depend upon the choice 
of the representation of the fields. 

1. Waveguide mode operators 

a. Mirror symmetry 

We start with the waveguide mode operators. The following general considerations apply 
to any single-mode waveguide that obeys mirror symmetry. For such a waveguide, one can 
always choose a mirror plane perpendicular to the waveguide such that the dielectric function 
e(x, y, z) satisfies 

^{-x^y^z) = e{x,y,z), (Al) 

where the x-axis is along the direction of the waveguide. From the Maxwell's equations 

d ^ 

V X E{x, y, z, t) = -/i— if (x, y, z, t), 

d ^ 

V X H{x, y, z, t) = +e(x, y, z)-q^E{x, y, z, t), (A2) 
it is straightforward to show that if the fields of the following form 

En{^,y,z,t) ^ le^'^^-— * ^ {A,{y, z), Ay{y, z), A^{y, z))e''^''-'-' 

= Eii{x,y,z)e-''^\ 
HR{x,y,z,t) = Ee*'^^— * = {B^{y, z), By{y, z), B,{y, z))e''^^--^^' 

= HR{x,y,z)e-'^' (A3) 
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is a solution of the Maxwell's equations for the waveguide, then the following fields are also 
a solution: 

Ei(x,2/,z,t) = {-A,{y,z),Ay{y,z),A,{y,z))e-'''^^-''^' 

HL{x,y,z,t) = {B,{y,z),-By{y,z),-B,{y,z))e-"'^^-"^' 

= HL{x,y,z)e-'^K (A4) 

The two sets of eigenmodes defined by Eqs. flA3l) and flA4p obey the mirror-symmetry 
transformation: 



(A5) 



where r = (x, y , 2;) , i.e., 



7r^{x, y,z) = {-x,y,z), 
T^x{Ex, Ey, Ez) = {—Ex,Ey,Ez), 

7T,{H,, Hy, Hz) = {H,, -Hy, -Hz) (A6) 

In the second-quantization form, one thus defines a mirror-transormation operator ir^ 
such that y 

(A7) 



n^c^j^{x)7r^ = c{{-x) 



Note that the mirror-transformation operator vr^ is both unitary and Hermitian 

-1 - _ ^ (A8) 



due to the choice of the relative phase of E^ with respect to Er in Eq. (lA4p |23j. Also, in 
the second quantization form, tt^ operates only on operators, not on numbers or functions. 

As an example, one can show that the mirror transformation of a right-moving photon is 
a left-moving photon: 

71^ / rfx/(x)c5j(x)|0) 



dxf{x) 

dxf{x)c[{-xm 
dxf{-x)c{{xm, 



7r.|0) 



(A9) 
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where, without losing generahty, the ground state is assumed to have positive parity: 7r^.|0) = 

10). 



b. Time-reversal symmetry 

From the Maxwell's equations Eqs. (lA2p . one can show that if T]Eji{r) and riHfi{r) is 
an eigenmode, where Eji and Hji are the fields in Eq. (]A3P and rj is an arbitrary complex 
number, then 



,* A* „—ikxX 



-r/*#* (r) = rj* i-B' ) e 



(AlO) 



is also an eigenmode 



24| . E*^{f) and —H'^if) are the time-reversed fields corresponding to 



Em 



Ejiif) and Hjiif). Since the waveguide is single- moded, we must have 

^L(r), 
-if^(f)=^z.(r), 

where we have chosen the proportionality constant to be 1. 
Thus we define a time-reversal operator T 

T [t^Ek] = 7]* El, 



(All) 



(A12) 



In the second quantization form, one defines a time-reversal operator T such that 

T(r/4(x))T-i = r/*c[(x). (A13) 

As an example, for any right-moving one-photon state, the time-reversed state is 

T ! dx/(x)4(x)|0) 



(AM) 



dxr{x) Tc^^{x)T-^ T|0) 

dxr{x)c[{x)m, 



where, without losing generality, we have assumed T|0) = |0). 

As a side remark, we note that combining the mirror and time-reversal symmetries allow 
us to make further statement on the fields. For example, can be chosen to be purely 
imaginary, while Ay and A^ can be chosen to be purely real. 
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2. Single- mode cavity 



For a single- mode cavity that obeys mirror symmetry, one must have 25 1 



±E{r), 



(A15) 



where the fields that transform with the upper(lower) sign is called an even(odd) mode. 
Thus, if a is the creation operator of the single-mode, one has 



(A16) 



where +(— ) sign is for even(odd) mode. Similarly, since a single mode must be mapped into 



itself, one has 



(A17) 



As an example, we now show that the symmetry considerations allow us to use a single 
real number to characterize the coupling between a single-mode cavity and the waveguide. 
The interaction between them can be written as 

Hi = j dx5{x) Vic^jf^^x)a + V*a)cR{x) + j dx5{x) V2c\{x)a + a) cl{x) . (A18) 

By writing down the interaction, we only require it to be Hermitian. 
The mirror-symmetry invar iance requires tTxHitTx = Hj: 



dx5{x) Vi'KxC^j^{x)-KxT^xa'Kx + V^i*vra;aV^7r^CR(x)7ri 



dx5{x) 



dx5{x){±) Vic\{-x)a + V^a^CL{-x) + dx6{x){±) V2c'j^{-x)a + V^a^cni-x) 



dx5{x) (±) \Vic{{x)a + V^*a^CL{x)] + / dx5{x) (±) \v2C^ji{x)a + V^a^cnix)] . (A19) 



Comparing with Eq. (lAlSp . one thus has ±Vi = V2. 
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The time- reversal invariance requires THjT ^ = Hj: 



THjT~^ 

dx6{x) \v{Tc^j^{x)T-^TaT-^ + ViTa^T~^TcR{x)T-^ 
+ j dx5{x) V^Tc{{x)T-^TaT-^ + V2Ta^T-^TcL{x)T~^ 

dx5{x) Vi <^i^{x)a + Vio^ cl{x) + / dx6{x) V2*c^(x)a + V2a^CR(x) 



(A20) 



Comparing with Eq. flAlSI) . one concludes that = V2. 

Combine the above results, one reaches the result that Vi = V2 = V{ = V for even 
mode, while Vi = — V2 = —V* = iV for odd mode, where ^ is a real number. Thus the 
coupling between a single-mode cavity and the waveguide could be represented by a single 
real number V, when the structure obeys mirror and time-reversal symmetry. 



3. WGM-type resonator 

We now consider a resonator which possesses azimuthal rotational symmetry and supports 
a pair of degenerate counter-propagating WGMs. Many resonators are of this type: two- 
dimensional dielectric cylinder, microsphere, microdisk, microtoroid, and ring resonator. 



a. Mirror symmetry 

We consider a two-dimensional cylindrically symmetric system. In the cylindrical coor- 
dinates (p, (f), z), if the fields of the following form 

Em{p,<P',z,t) ^ ie*™^'— * ^ {A,{p,z),A^{p,z),A{p,z))e'"^'^'-'^' 
= E^{p,ct>\z)e-'^\ 

= Hm{pA',z)e-"'' (A21) 
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are a set of solutions of Maxwell's equations for the resonator, then the following fields are 
also a set of solutions : 

H_UP,<P',z,t) = {-B,{p, z),B^{p, z),-B,{p,z))e-'"''^'-^' 

= H^^{p,(t>\z)e-'^K (A22) 

The unit vectors in the cylindrical coordinates could be written in the Cartesian coordi- 
nates as 

p = cos 0' X + sin 0' y, 



sin (/)' X + cos (j)' y\ 



(A23) 



where the a:'-axis is the azimuthal origin of 0' = 0, and needs not to be coincident with 
X-axis that is along the direction of the waveguide. 

Thus, the two sets of eigenmodes of Eq. flA21l) and Eq. flA22p in the Cartesian coordinates 
could be expressed as 

Era{r) = [{Ap COS (j)' - A^ sin 0') x' + {Ap sin (p' + A^ cos </>') y' + AJ] e'"^^' , 

Hm{r) = [{Bp cos 0' - B^ sin 0') x' + {Bp sin 0' + B^ cos 0') y' + BJ] e'""^' , (A24) 

and 

E^rn{r) = [{Ap COS (p' + A^sin^') x' + {Ap sin (p' - A^coscp') y' + AJ] e''""^' , 

H_^{r) = [{-BpCos(j)' - B^smcj)') x' + {-BpSmcj)' + 5^cos0') y' - BJ] e"*™"^', (A25) 

respectively. 

An azimuthal rotational symmetric structure also possesses mirror symmetry. The choice 
of mirror plane relative to the x'-axis however affects the specific form of transformation 
between the fields. Here we examine this degree of freedom. 

As before, let x-axis be along the waveguide direction, and the mirror plane is the y-z 
plane that perpendicular to the x-axis. If the angle between the x'-axis and the x-axis is (po 
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(see Fig. [T71) . the bases of Eqs. (1A23P in x — y coordinate become 



p = cos{(f)' + (f)o)x + sm{(j)' + y, 



sin(0' + (f)Q)x + cos(0' + 0o) y, 



z = z. 



(A26) 



and the fields become 



im(j) —im(f>(] 



Em = [{Ap COS 0-^0 sin 0) x + {Ap sin (f) + A^ cos (f)) y + A^z\ e*""^e 
Em = [{Bp cos (p-B^ sin 0) £ + (5^ sin (j) + B^ cos 0) y + e^'^'^e-*'"'^^, 



[{Ap cos + sin (f)) x + {Ap sin - cos (f)) y + A^z] e 



= [{-Bp cos 0-5,^ sin 0) x + (-5^ sin + 5,^ cos 0) y - 5^1] e-*'"<?'e+*'"'^" , (A27) 

where = 0' + 0o is the azimuthal angle of the point r in the x-y coordinate. The mirror 
transformation that maps x — > —x is equivalent to the mapping ^ tt — 0, while with p 
and z unchanged. The fields thus transform according to 



Em (vT^r) 



Hm (vT^r) 



.l^rn^-2im^0ff_^ (f) , 



(A28) 



where the transformation rules of on fields are given by Eq. flA6 
In the second quantization form, accordingly one has 



(A29) 



where = 7r/2 — 0o is the angle between the x'-axis and the mirror plane {y-z plane). In 
the special case where the x'-axis and the x-axis are coincident (0o = 0), one has 



IT T a' IT. 



-irb\ 



TT. 



(A30) 



where in the special case where the x'-axis and the y-axis are coincident (0o = 7r/2), one has 



(A31) 
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b. Time-reversal symmetry 



The discussion of the transformation by the time-reversal operator T is the similar to 
that of the waveguide case. If the fields riE„^{f) and rjHmir) are a set of solutions, then 
i]*E^{r) and —ri*H^{f) are also a set of solutions. From the functional form, one concludes 
that 

-H*Jr) = H.^{r), (A32) 

where we have chosen the proportionality constant to be 1. We thus define a time- reversal 
operator T such that 

T r]Em = r]*E-ni, 

T [riHm] = V*H.m. (A33) 

In the second quantization form, one defines a time-reversal operator such that 

T {7]a^) = r]*bl (A34) 

As a side remark, we note that combining the mirror and time-reversal symmetries allow 
us to make further statement on the fields. For example, can be chosen to be purely 
imaginary, while Ap, and A^ can be chosen to be purely real. 

c. Waveguide-ring resonator interactions 

The general form of the waveguide-ring resonator interaction can be written as 

Hj= [ dx6{x)[VaC^j^{x)a + V*a^CR{x)]+ [ dx 6{x)[Vbc{{x)b + V^*b^CL{x)], (A35) 



where we only require it to be Hermitian. The symmetry invariance further impose con- 
straints on the coupling constants: 

The mirror-symmetry invariance requires Ti^Hjir^ = Hj: 

= j dx5{x)[Vacl{-x)e-''''^h + V:e+''''^h'^CL{-x)] 

+ /c/x5(x)[H4(-x)e+'2'"^a + V;e-'2'"Vcij(-x)], (A36) 
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which yields VbC*^™''^ = Va, where (f) = 7r/2 — (f)Q, as previously defined. 
The time- reversal symmetry invariance requires THjT"^ = Hf. 

THiT~^ 

= j dx5{x)[V:c{{x)h + Vah^CL{x)] + j dx5{x)[VCc\{x)a + V^a)cR{x)l (A37) 

which yields Va = V^ . 

Combine the results together, one has 

Va = Ve'"^^, 

V = Ve-'""^, (A38) 

where V = \Va\ = \Vb\ is a real number. 



APPENDIX B: RESONATOR- ATOM INTERACTIONS 

In this section we outline the procedures to obtain the interaction between the resonator 
and the atom in the second quantized form. This derivation gives the general form of the 
interaction. 

The interaction between the resonator and the atom is given by 

Hi = -d- E, (Bl) 

where d = er is the electric dipole moment of the atom (e > is the magnitude of charge), 
and E is the electric field of the resonator modes. This expression is valid when the wave- 
length of interest is much longer than the dimension of the atom. The field is evaluated at 
the location of the dipole r = tq. 

The atomic electric dipole operator could be expanded using the bases of the atom: 

d = (|-)(-| + |+)(+|)d(|-)(-| + |+)(+|) 
= (l+)(+|d|-)(-|) + (|-)(-|d|+)(+|) 
= d\+){-\+d*\-){+\ 

= da+ + (f CT_, (B2) 

— * 

with d = e(+|r| — ). The terms (+|d|+) = (— |d| — ) vanish since d has odd parity. Here the 
atomic basis |— ) and |+) are taken as eigenstates of parity. Thus the atomic electric dipole 
operator is entirely off-diagonal in the atomic bases. 
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The electric field operator is expanded using the eigenmodes of the resonator. We will 
treat two cases here: the single-mode cavity and a ring-resonator. 
For a single-mode cavity, the electric field operator is 

E(r) = a(^(r) + aV*(r), (B3) 

where is the field profile of the mode, up to a global phase that does not affect the physics. 
Consequently, the interaction takes the following form: 

Hi = - (^a(f{r) + aV*(r)) ■ + dV„) 

= - (^0(r) ■ (tj aa+ - (^0*(r) ■ cf ) - (^0(r) ■ cf ) aa_ - (^0*(r) ■ dj aV_ 

= gaa+ + g*a^a-, (B4) 

where g = — 0(r) ■ d. We have omitted the latter two terms which give zero contribution 
for single-photon transport. For single-photon transport, the rotating wave approximation 
(RWA) becomes exact. 

For a resonator that support a pair of degenerate WGMs, the electric field operator is 
given by 

E(r) = aMr) + a^^ir) + 60,(r) + b^$l{r) 

= (a + 6t)0(r) + (at + 5)0*(r), (B5) 

where we have used 0a (r) = 0^(1") = 0(r)- Thus, the interaction is given by 



Hr 



{a + b^) 0(r) + (a^ + 6)0*( 



= -0 ■ d{a + 6^) (T+ - ■ d* [a + b^) a_ - 0* ■ d{a^ + 6) a+ - 0* ■ d* {a^ + b) 
= —0 ■ daa^ — (f) ■ d* b'^a^ — 0* ■ dba^ — (p* ■ d* a)a^ 

= {gaaa+ + gla^aJ) + {ghba+ + glb^aJ) , (B6) 

where ga = —4>-d, and gh = —0* ■ d. This gives the most general form of the interaction 
for a WGM-type resonator interacting with a two-level atom. Note that in systems without 
underlying time-reversal symmetry, \ga\ 7^ \gb\- 

Here we make further remarks on the form of the interaction, Eq. (\B6\i : 

1. If the transition from state |— ) to |-|-) corresponds to a Am = of a real atom (m is 
the magnetic quantum number), ci is a real vector. On the other hand, for an atomic 
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Am = ±1 transition, such as might be induced by circularly polarized light, d is 



271, 



28| . When ci is a real vector, one has ga = gl = g, 



necessarily a complex vector 26 
and Eq. ( ]B6[) becomes 

Hj = {gaa+ + g*a^ciJ) + {g*ha+ + gh^aJ) . (B7) 

For the kind of experimental setup in Refs. , if the WGMs are TE mode so 

the electric field is linearly polarized along the z-axis (please refer to Fig. [T7] for axis 
orientation), such an electric field will induce only Am = transitions between atomic 
states. This yields a real dipole moment d and thus one has ga = 9l- 

2. Other simplifications of Eq. ( ]B6[) can be made. For example, using the fact that in 
the cylindrical coordinate system (p, 0, z) one has 



{Ep, -iE^, E^) e 



-im(f> 



{B8) 



(see Appendix [XI or Ref. (29|), and suppose that the electric dipole of the quantum 
dot is transverse to the (j) direction, and the azimuthal origin lying at the location of 
the quantum dot, the interaction takes the form 



Hj = [gaa^ + g*a^aJ) + [ghaj^. + g*h'^aJ) 
where ga = 9b = 9- Such a choice was made in Ref. 6| 



(B9) 



APPENDIX C: WAVEGUIDE-RING RESONATOR SYSTEM 

Before coupling to the atom, one has to calibrate the system of waveguide-resonator, and 
to quantify the physical parameters such as the coupling strength, the intrinsic loss, and the 
Q-factor of the resonator. Moreover, the waveguide-ring resonator system itself has been 
widely adapted in opto-electronic devices to achieves functions such as add-drop filter, delay 
line, and signal switching. In this appendix we give the transmission amplitude for this 
important special case. 

The transmission amplitude of a system of waveguide-ring resonator is given by 

m = ) 72 ' (CI) 
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which can be easily obtained by either starting with the Hamiltonian H in Eq. ([2]) and the 
state le"*") in Eq. ([7]) by omitting the atomic parts; or from the expression of t in Eq. fllObp 



by letting the atomic transition frequency Vt oo and the resonator-atom coupling — 0, 
which effectively decouples the atom from the waveguide-resonator part. Note only the 
magnitude of h appears in the amplitude, its phase does not. For simplicity, we assumed 
Va = Vb = V, and ga = gt = g- 

For the special case when h = 0, Eq. (1C1|) and other amplitudes yield the well-known 
simple forms: 

UJ — ujf. + i— — iT 

UJ — Uc + t— + ti. 

eaiu;) = (C2b) 

r{uj) = = eb{uj). (C2c) 

Note that the reflection amplitude r is always zero even when the transmission amplitude 
is not of magnitude 1. Only the phase-matched resonator WGM (a) is excited. It is well- 
known that the ring resonator acts as an all-pass filter when the resonator's intrinsic loss 
(I/tc) is negligible such that I/tc ^ F. In this limit, \t\ ^ 1 for all frequency uj. On the 
other hand, when the resonator intrinsic loss dominates such that l/r^ ^ F, one again has 
|t| ~ 1 for all frequency u, and the ring resonator again acts as an all-pass filter. 

On the other hand, for the lossless case I/tc = 0, the locations of the local minima and 
maxima of the transmission spectrum can be obtained by solving for u from = 0, which 
yields 

UJ = ujc, Uc ± \J \hf — r"^ = u±. (C3) 
The number of dips in transmission thus depends upon the sign of 

\hf-T\ (C4) 

i.e., if \h\ > F the spectrum exhibits a doublet structure; while if \h\ < F, the spectrum 
shows a single dip. 
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APPENDIX D: REFLECTION AND COUTER-PROPAGATING WGM EXCI- 
TATION SPECTRUM 

Here we plot the reflection spectrum R = \r{uj)f and the counter-propagating WGM 
excitation spectrum \ei,f in Fig. [T51 and Fig. [111 The parameters are the same as those in 
the transmission spectrum plot in Fig. |H 



APPENDIX E: WAVEGUIDE-RING RESONATOR- ATOM SYSTEM AT h = 



In this section, we give the results wherein the inter-mode backscattering is not present. 
The results are written in the form revealing the underlying physics. 
When h = 0, solving Eqs. (!9al) - (!9el) for Q, t, r, ea, e^, and Cq gives: 



Q 



(^uj-n + (u-Lj, + i^^ - 2\g\ 



iT 



uj — Q 



219 



2 uj — uJc + i— + iT' 



{to -^1 + i^)iuj -uJc + i^ + iT) - 2\g\ 



+ 



iT 



UJ — Uc + i— + iT 



ea 



eb 



+ 



2 



(u — Q + i—)(LV — uJc + i— + iT) — 2|fifp uj — Uc + i— + iT 



2 



(u; - n + i^){yj - + + iT) - 2\g\^ 
"gV* 



u ~ Uc + i — h iT 



(Ela) 
(Elb) 

(Elc) 

(Eld) 

(Ele) 
(Elf) 



" {uj -VL + i^){uo -uJc + i^^+iT) -2\g\^' 

which are valid in both strong and weak coupling regimes, g = \g\e^^^. Note that although 
t, Ca, and Cq depend upon the magnitude of g only, both r and e^, however, are proportional 
to g^ and thus depend upon the phase of (7 as ~ e^*^". Only when h = 0, can one eliminate 
the phase dependence on g by properly choosing the azimuthal origin of coordinate system. 
The above expressions give explicitly the phase dependence without making the assumption 
that g is real. 

The form of the amplitudes, Eqs. (lElbp - (lElf[) . is a direct consequence of the fact that 
the ring resonator supports two degenerate whispering-gallery modes. Two orthogonal non- 
propagating modes can be formed as linear superpositions of these two degenerate modes: 
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one of them has non-zero amphtude at the location of the atom, and gives a contribution 



which resembles that of the waveguide-single-mode cavity-atom system [30|, but with an ef- 
fective atom-resonator coupling \plg and an effective resonator- waveguide coupling |V"|/-\/2. 
The other mode on the other hand has zero amplitude at the location of the atom, and 
thus decouples from the atom. This mode gives a contribution which resembles that of 
the waveguide-single- mode cavity system, with an effective resonator-waveguide coupling 
|\^|/v^- This explains the form of the amplitudes. 
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FIG. 1: (Color online) Schematics of the system. The single-mode waveguide is denoted by the 
blue channel. The ring resonator is denoted by the green ring. The black dot denotes the two- 
level atom, of which the two energy levels are plotted on the right. The waveguiding modes are 
described by c^^{x) and cj^(x). The two degenerate whispering gallery modes are described by 
(counter-clockwise) and 6^ (clockwise). 
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FIG. 2: (Color online) Fitting to experimental data of system of waveguide-microtoroidal resonator, 
r = \V\^/2vg. (a) Lower curve: |/i|/27r = 5.50092 MHz, 1/tc/2tt = 7.49983 MHz. r/27r = 6.89887 



MHz. F'^ = 1.34818'^ Y -|- (I/tc)^, which is somewhat close to the critical coupling condition. 
Upper curve: under-coupling. Blue curve: \h\/2TT = 7.57069 MHz, 1/tc/2tt = 8.4642 MHz. F/27r = 
1.13735 MHz. F^ < + (I/tc)^. The overall decay rate k in Ref. [5| is equal to T + l/r^ in 
this article, (b) |/i|/27r = 7.64947 MHz, l/rc/27r = 0.250879 MHz, F/27r = 0.0194301 MHz. 
F^ <^ -|- (I/tc)^ and \h\ » 1/tc S> F. The numerical values are obtained by least squares 
method using the amplitude of Eq. ([13]) . 
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FIG. 3: (Color online) The transmission spectrum T (blue curve), reflection spectrum R (red 
curve), and the atom excitation spectrum A = |eq|^ (green curve) of the coupled waveguide- 
microdisk-quantum dot system. The normalized spectrum 2^|eg|^ is plotted. Upper panel ((al)- 
(cl)): dephasing 7p = 0. Lower panel ((a2)-(c2)): 7p / 0. (al), (a2): VL = uj^. (bl), (b2): 
Vl — uj(. = \h\. (cl), (c2): fi — = —\h\. Parameters used: g/2TT = 6, h/2'K = —9.6, l/Tg/27r = 0.16, 
1/tc/2tt = 0.76, r/27r = 0.44, and 7p/27r = 2.4. 
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Transmission spectrunn T 




FIG. 4: Map of the transmission spectrum T as g and h are varied, h is real for this map. Also, 
9a = 9b and g = \ga\- The location oi lo = u>c is indicated by the gray tick in each figure. 
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FIG. 5: Map of the group delay d<j)/du> as g and h are varied, h is real for this map. Also, ga = gt 
and g = \ga\- The location uj = uJc is indicated by the gray tick in each figure. 
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FIG. 6: Map of the atomic excitation spectrum \eq\ as g and h are varied, h is real for this map. 

r I 1 2 

Also, ga = gt and g = \ga\- The normalized spectrum ^jrl^?! plotted. The location of uj = Uc is 
indicated by the gray tick in each figure. 
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Phase matched WGM exciation 
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FIG. 7: Map of the phase- matched WGM excitation \ea\ g and h are varied, h is real for this 
map. Also, ga = gb and g = \ga\- The normalized spectrum 2^|ea|^ is plotted. The location of 
u; = is indicated by the gray tick in each figure. 
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Anti-crossing of resonances 
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FIG. 8: Anti-crossing between the atomic and cavity resonances as h is continuously varied, g = 
2.5T. This depicts the evolution of the resonances in the column g = 2.5T in Fig. [H 



Transmission spectrum with intrinsic losses 
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FIG. 9: (Color online) Map of the transmission spectrum with intrinsic losses as g and h are 
varied. Red curves: only atom dissipation is non-zero (l/rg = T and 1/tc = 0). Blue curves: 
only resonator dissipation is non-zero (1/tc = T and l/r^ = 0). The intrinsic atom dissipation 
only affects resonances with atomic nature, and the intrinsic resonator dissipations only affects 
resonances with cavity nature, h is real for this map. Also, ga = gb and g = \ga\- The location of 
a; = is indicated by the gray tick in each figure. 
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Transmission spectra with intrinsic losses (g = 2.5, h = 10) 
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FIG. 10: (Color online) The transmission spectrum for fixed g and h as the intrinsic atom and 
resonator losses are increased, respectively. From left to right: 0, O.SF, F, 5F, lOF, and 20F. 
Resonances of large cavity nature are less affected by the atom dissipation; similarly, resonances of 
large atomic nature are less affected by the cavity dissipation. The location of w = is indicated 
by the gray tick in each figure. 
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Transmission spectra (g = 2.5, h = 2.5e* ^0 
1— r— 1 — X r — 1— ^ / — 1 








photon angular frequency 



FIG. 11: The transmission spectrum when h is complex, g = 2.5T. h = 2.5e^^'^T. From left to 
right: 6 = 0, 7r/4, ir/2, 37r/4, and vr. The spectrum of = is mirror image to that of = vr, and 
6 = 7r/4 to Sit/ 4. The spectrum of ^ = 7r/2 is symmetric with respect to u; = tOc- The location of 
LO = LOc is indicated by the gray tick in each figure. 
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FIG. 12: The general symmetry property of the transmission spectrum when A 7^ 0. Left: Q — u>c = 
2r, h = he'^'I'^V. Right: - = -2r, h = Se^^'^/^F. l/r^ = 2F, l/r^ = F, and g = 3T for both 
plots. 
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Transmission spectrum T for detuned case 

\^ 0.3 2.5 a 

A V 




FIG. 13: Transmission spectrum for the lossless detuned case, h = 0, ga = gt and g = \ga\. The 
location of uj = u>c is indicated by the gray tick in each figure. 
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Transmission spectrum witfi intrinsic losses 




FIG. 14: (Color online) Transmission spectrum for the detuned case with losses. Red curves: only 
atom dissipation is non-zero (I/t^ = F and 1/tc = 0). Blue curves: only resonator dissipation is 
non-zero (1/tc = F and l/r^ = 0). h = 0, ga = Qb and g = \ga\. The location oi u = ujc is indicated 
by the gray tick in each figure. 
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Reflection spectrum R 




FIG. 15: Map of the reflection spectrum R as g and h are varied, h is real for this map. Also, 
5a = gb and g = \ga\- R + T = 1. The location of lo = uJc is indicated by the gray tick in each 
figure. 
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FIG. 16: Map of the counter-propagating WGM excitation \eb\ as g and h are varied, /i is real for 
this map. Also, Qa = gt and g = \ga\- The normalized spectrum 2^|ef,|^ is plotted. The location 
to = LUc is indicated by the gray tick in each figure. 
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FIG. 17: (Color online) The orientation between the x'-y' and x-y coordinate systems. The angle 
between the x'-axis and x-axis is <^q. The angle between the x'-axis and the mirror-plane 
plane) is thus (j) = n /2 — 
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